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In the classical Josephson effect the phase difference across the junction is well 
defined, and the supercurrent is reduced only weakly by phase diffusion. For meso- 
scopic junctions with small capacitance the phase undergoes large quantum fluctu- 
ations, and the current is also decreased by Coulomb blockade effects. We discuss 
the behavior of the current-voltage characteristics in a large range of parameters 
comprising the phase diffusion regime with coherent Josephson current as well as 
the supercurrent peak due to incoherent Cooper pair tunneling in the Coulomb 
blockade regime. 



1. Introduction 

The Josephson effect beautifully embodies quantum tunneling and superconductivity (I], || . Joseph- 
son systems have been a subject of fundamental research for decades, which have led to useful 
applications, e.g., in highly sensitive magnetometry and metrology. Like in the Ginzburg-Landau 
theory for bulk superconductivity, where the phase of the order parameter is considered as a qua- 
siclassical variable, the phase difference across a Josephson junctions can frequently be treated 
classically. However, in the last ten years or so, new lithography and low-temperature techniques 
have allowed the fabrication and measurement of small Josephson junctions affected by the capac- 
itive charging energy of single Cooper pairs ||, ||] . A large charging energy will render the Cooper 
pair number on either side of the junction classical and thus cause large quantum fluctuations of 
the conjugate phase variable. 

In this article we show how the classical Josephson effect, that is a supercurrent / flowing at 
vanishing voltage V, gradually evolves via the classical phase diffusion regime into the supercurrent 
peak in the Coulomb blockade regime where Cooper pairs tunnel incoherently across the Josephson 
junction. After a brief review of the classical Josephson effect in sec. || and phase diffusion in sec. || 
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Fig. 1. (a) Schematic view of a Josephson junction where two superconductors are separated by an insulating 
barrier, (b) Equivalent circuit with Josephson coupling Ej and capacitance C. The voltage across the junction is 
denoted by Vj. 



we discuss Cooper pair tunneling in the Coulomb blockade regime in sec. [|. The turnover between 
the classical and the quantum regimes is discussed in sec. [| Finally, sec. ^ contains some concluding 
remarks. 



2. Classical Josephson Effect 

A Josephson junction consists of two superconductors separated by a thin insulating barrier as 
shown in fig. [l]a. Its dynamics may be described in terms of the phase difference ip of the two 
condensate wave functions on the left and right sides of the junction and the number n of Cooper 
pairs on the capacitor formed by the junction. These operators n and (p obey the commutation 
relation [n, ip] — — i, and the junction Hamiltonian may be expressed as 

Hj = E c n 2 - Ej cos(<p). (1) 

The first term describes the charging energy associated with the capacitance C of the tunnel junction 
where E c — 2e 2 /C is the charging energy corresponding to a single Cooper pair. The second term 
arises from the tunneling of Cooper pairs through the junction characterized by the Josephson 
coupling energy Ej. These properties may be expressed in terms of the circuit shown in fig. 
The Josephson relations link the phase to the voltage across the junction 

V, = & (2) 

and to the supercurrent 

I = I c sm(p). (3) 

The critical current I c is related to the Josephson coupling energy by I c = 2eEj/h. These relations 
allow for a supercurrent flowing at zero voltage if ip remains constant in time. 

3. Phase Diffusion in Josephson Junctions 

In the real world Josephson junctions are coupled to an electromagnetic environment which may 
be described by an impedance Z(oj). For simplicity, we will in the following mostly consider the case 



Mesoscopic Josephson Effect 




Fig. 2. Circuit modelling a real experimental set-up: a Josephson junction characterized by the Josephson energy 
Ej and capacitance C is coupled to a voltage source V via a resistor R. The voltage drops across the junction and 
the resistor are Vj and Vr, respectively. 



of an ohmic resistor [Z(lj) = R] but a more general case will be addressed in sec. [|. The resistor 
will give rise to Nyquist noise and therefore to a diffusive behavior of the phase difference <p. 

To be specific, we consider the circuit shown in fig. || where the Josephson junction is coupled 
to the voltage source via an ohmic resistor. The voltages across the junction and the resistor are 
denoted be Vj and Vr, respectively. It is useful to introduce the dimensionless resistance p = R/Rq 
where Rq = h/4e 2 is the resistance quantum for Cooper pairs. 

In the classical limit this system can be described by the Langevin equation jl], ^| 

c (SK(S)+'*x- <«) 

The term on the righthandside represents the Nyquist noise with (<$Vr,(£ + t)<5Vr (£)} = (2R/(3)8(t) 
for a resistor at inverse temperature = 1/k-sT. Here, SVr denotes the fluctuation of the voltage 
about its average value. 

This problem has been solved by Ivanchenko and Zil'berman || in the overdamped limit where 
2ir 2 p 2 Ej -C E c . The stationary solution of the Fokker-Planck equation corresponding to eq. (|4|) can 
be expressed in terms of modified Bessel functions of complex order and yields the Cooper pair 
current || 



Ic \ I-iPE 3 v{pEj) 

as a function of the dimensionless applied voltage v = V/RI C . While the supercurrent at zero voltage 
is destroyed by the fluctuations causing the phase diffusion, the current-voltage characteristics (||) 
display peaks at small but finite voltages as shown in fig. ||a. With decreasing temperature the I-V- 
curve (ra) becomes closer to an ohmic line for voltages up to RI C , while the peak in the I- Vj -curve 
in fig. 0b moves towards zero voltage. 

We analyze the result (p|) by considering the zero bias differential resistance 



R °-~dT 



dV 
~dl 



- R. (6) 

v=o 



Vj=0 

Ro is defined with respect to the junction voltage Vj, and 1/Ro describes the slope at the origin in 
fig. From eq. (||) one obtains 

^ = I (7) 
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For /3Ej 3> 1, the I-V-curve is very close to / = V/R with an exponentially small difference 

^=27r/J£?jexp(-2/3Ej). (8) 

This shows that for temperatures much lower than the height of the periodic potential for (p in 
cq. ([!]), the dynamics of the phase is thermally activated. 

On the other hand, for (3Ej <C 1 one finds in the overdamped limit 

^ = 2 (9) 

In this case, the I-V-curve becomes 

I 2 RV 

2 V 2 + {2eR/hp) 2 { ' 

which corresponds to a broad peak structure. 



4. Coulomb Blockade in Josephson Junctions 



We now turn to the case of ultrasmall tunnel junctions with spatial dimensions so small that 
E c ^> Ej, Then the phase can no longer be treated as a quasiclassical variable, rather the charge 
on the junction capacitance will approximately follow classical statistics. As we have seen above, 
due to phase fluctuations, a supercurrcnt at zero voltage is no longer possible. However, a Cooper 
pair current may flow at finite voltages if the electromagnetic environment is able to absorb the 
energy 2eV gained by a Cooper pair tunneling through the junction. In contrast to sec. |^, it is now 
necessary the introduce the external impedance on the quantum level. Due to its linearity the exter- 
nal impedance may be modeled by a possibly infinite number of LC-oscillators. The corresponding 
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Hamiltonian then reads 



-ffimp — 



2C n \2ey 2L ra 



-(VR ~ Vnf 



The external impedance in terms of the inductances L n and capacitances C n is given by 

-i -l 



Z(w) 



/* _ _ 

/ dte-^ V 

A n=l 



cos(o> n i) 



(11) 



(12) 



where w„ = (L n C n ) 1//2 . The coupling between the external impedance and the phase difference ip 
appears through the phase 

<PR=*J<&{V-V 3 ) = ^Vt-<p. (13) 
This phase may formally be attributed to the resistance by making use of the Josephson relation 



Calculating the Cooper pair current from the Hamiltonian 

H = Hj + -ffimp 

perturbatively to lowest order in Ej, one obtains |l| 



I = 



[P(2eV) 



P(-2eV)] 



(14) 



(15) 



Here, P(2eV) and P(—2eV) are the probabilities that the energy 2eV of the tunneling Cooper pair 
is absorbed or provided by the environment, respectively. This probability depends on the external 
impedance and on temperature through B 



with 



2nn 



J(t) = 2 



dt exp 



J(t) + T Et 
n 



dui ReZ t (oj) e 



1 



CO 



Rq 1 

This expression contains the total impedance seen by the Josephson junction 

Zt(u>) ' 



iluC + Z-\io) 



(16) 



(17) 



(18) 



which is given by the external impedance in parallel with the capacitance. 

For an ohmic environment [Z(u>) = R = pRq] the result ( |l5| ) may be evaluated further. At zero 
temperature one finds a zero bias anomaly given by I ~ V 2p ~ x . This shows that at T = the 
perturbative result (|l^) is only valid for large p > 1/2 where the current remains small for all 
voltages. In this case the I-V-cmve displays a peak near voltages of order e/C as shown in fig. [l| 
This can easily be understood from a simple Coulomb blockade picture. At V — e/C the voltage 
gain 2eV just equals the charging energy E c . 

However, the experimentally relevant situation is an environment of low impedance (p -C 1) pOf . 
On the one hand, it is difficult to fabricate large resistances at frequencies around E c /%, which 
typically is of the order of 10GHz, and on the other hand the resistance should not be too large 
to avoid heating. Furthermore, since charging effects are washed out by thermal fluctuations, we 
are interested in the regime of low but finite temperatures (3E C 3> 1. To proceed for parameters in 
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CV/e 



Fig. 4. Cooper pair current-voltage characteristics I(V) in the Coulomb blockade regime at zero temperature. The 
peak is sharpening with increasing p = 2, 20, 100. 



this range, we first note that the function J{t) is equivalent to the position autocorrelation function 
of a free damped particle and can be evaluated in closed form for ohmic damping jll). For long 
times, J {t) describes diffusive behavior linear in time for finite temperatures and logarithmic in 
time at zero temperature. We therefore neglect exponentially decaying terms. This restricts us to 
low voltages which is the regime of interest. Assuming furthermore p <C (3E C in accordance with the 
above considerations, we arrive at [Q 



J(t) 



-2p \ln 



(3E C 
ir z p 



sinh 



nt 



+ 7 + i^-sign(t) 



(19) 



Here, 7 = 0.5772 ... is the Euler constant. Inserting ( |l9| ) into ( |l6| ) and making use of (|l5| ) one finally 
obtains H || @ 



E, 



V 2tt 2 



tt Ej 2p (pE c V- 2p 



|r(p-i 



r(2p) 



sinh(/3eF). 



(20) 



This perturbative result is only correct if the current is not too large which may be the case for 
very small temperatures and very small damping. In fact, as noted above, at zero temperature one 
finds a divergence at low voltages. Therefore, the validity is restricted to not too low temperatures 
/3Ej <C p. Since E c » Ej, this together with the upper bound on temperature \(3E C » p] still leaves 
a rather large range where (|2(]) is applicable. 

The current-voltage characteristics ( p0| ) is shown in fig. ^ (see sec. ||) as dashed line for f3Ej = 
2, p = 0.04 and (3E C = 1 and 10 5 . It exhibits a peak at finite voltage which for small p is given by 



K ia x = ^(l+4C(3)p 3 + ...). 

ep 



(21) 



Here £(3) = 1.202... is a Riemann number. In the limit of small p, high temperatures and for 
(3eV <C 1 one recovers the result of classical phase diffusion (|l0|). 

A nonohmic environment of practical interest is a finite LC transmission line terminated by an 
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Fig. 5. Cooper pair current-voltage characteristic for an ultrasmall tunnel junction coupled to a finite transmission 
line with R^/Rq = 0.01, Rl/Roo = 0.1, and ujqRocC = 1 at zero temperature. The peaks correspond to excitations 
denoted by (N1N2N3N4), where N k is the number of quanta of the fc-th mode excited, a: (1000), b: (2000), c: (3000), 
d: (0100), e: (1100), f: (2100), g: (0010), h: (1010), i: (0001), j: (1001). 



ohmic load resistance Rl ]l5| , [16] . The transmission line is characterized by its resistance at infinite 
length Roc = (Lq/Cq) 1 ^ 2 where Lq and Co are specific inductance and capacitance per unit length. 
As a parameter describing the finite length I of the transmission line we choose the A/4- frequency 
wo = (n/2)u/£ where u = {LqCq)^ 1 / 2 is the velocity of wave propagation on the line. 

Depending on the ratio between load resistance and resistance of the infinite transmission line, 
the environment behaves quite differently. For Rl = Roo the impedance matching results in an 
ohmic impedance. On the other hand, for Rl very different from R^, the impedance displays sharp 
resonances. According to the above discussion, the properties of the external impedance should 
show up in the probability to absorb or emit energy i.e. in P(E), and therefore according to dl5| ) in 
the Cooper pair current-voltage characteristics. This is indeed the case as can be seen from fig. [| 
There we show the zero temperature Cooper pair current for Rl/Roo — 0.1 obtained from a direct 
numerical evaluation of eqs. (|l5|)-(p7|). The resonance peaks can be identified as single or multiple 
excitations of various transmission line modes ]l2fl . The measurement of the Cooper pair current 
thus allows for a spectroscopy of the environmental modes. 



5. Supercurrent Peak: Between the Phase Diffusion and Coulomb 

Blockade Regimes 

So far, we have discussed the behavior of a Josephson junction in the classical overdamped limit 
(sec. [3]) and in the quantum regime for charging energies much larger than the Josephson cou- 
pling energy (sec. ^). We now want to bridge the gap between these two limits by extending the 
perturbation theory presented in the previous section to infinite order in Ej. 

The summation of the perturbation series in Ej is only possible for an appropriate choice of the 
parameter regime. Motivated by the above discussion, we will continue to consider the overdamped 
regime where lor — 1/RC is much larger than the Josephson frequency wj = (2e/h)RI c . This is 
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identical to the condition 2ir 2 p 2 Ej <C E c given in sec. |^. In addition, the correlation function J(t) 
contains a thermal frequency scale given by v — 2tt /Ti(3 as can be seen from the last denominator 
in eq. (^). In the following we assume that wj -C v so that we may neglect terms exponentially 
decaying in time with v or faster. This results in the approximation 



J(t) = -2p 



JL|t| + 5 + i|si g n (i ) 



(22) 



valid for sufficiently high temperatures where pf3Ej <C 1. Further, we have introduced the abbrevi 
at ion 



S = 7 + ^ + V> 



I3E C ' T \2tt 2 P 



(23) 



where ip(x) denotes the logarithmic derivative of the gamma function. 

For the correlation function ( p2|) t he Cooper pair current can be evaluated exactly and expressed 
in terms of a continued fraction ]17| 



/ = 7 c Re 



sin(7rp) exp(— 2pS) 



2irp 



(24) 



with coefficients 



/ (3Ej \ 2 sm(npn) sin(irp(n + 1)) exp(— 2pS) 
\2np J n(n + l)(n - ivf3Ej)(n + 1 - ivf3E 3 ) ' 



(25) 



Since the continued fraction converges rapidly, we may linearize the sine functions appearing in 
cq. ( p5| ) for sufficiently small p. Then, the continued fraction may be evaluated with the help of a 
matrix recursion and one finds (l?]] 



2c 



E* r Im 



Il-if3eV/ivp{[3Ej) 
I-if3eV/ivp(PEj) 



with an effective Josephson energy 



E* 3 = Ej cxp(-p[^(l + ^ R /2vr) + 7 ]). 



(26) 



(27) 



Here we have rewritten the quantity S defined in eq. (23) in terms of the frequency lor and have 
neglected a term Trp/hflujR <JC 1 in the overdamped limit considered. Even though p has to be small 
for the expression ( |27| ) to hold, the correction to the bare Josephson energy may be substantial 
because for low temperatures the ip function grows logarithmically and may become large. On the 
other hand, in the classical limit where fi — > and e — > such that the flux quantum h/2e remains 
constant, the effective Josephson energy Ej coincides with the bare Josephson energy Ej. We thus 
recover in the classical overdamped limit the result (^|). 

In fig. ^| we show how the result (^6|) bridges between the phase diffusion result (^) shown as 



dotted line and the Coulomb blockade result (20) depicted as dashed line. The current-voltage 
characteristics have been calculated for f3Ej = 2 and p = 0.04. For large (3E C where charging 
effects should be important one obtains very good agreement with the Coulomb blockade result. 
With decreasing (3E C the current-voltage characteristics evolve differently and a crossover to the 
Ivanchcnko-ZiPberman result is found. 
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Fig. 6. The current-voltage characteristic of a Josephson junction with Josephson energy /3Ej = 2 is shown for 
charging energies /3E C = 1 and 10 5 and external resistance a = 0.04. The full line corresponds to our result (E4J) w hile 
the dotted line gives the standard Ivanchenko-Zil'berman (0) result and the dashed line depicts the prediction (201) 
for Coulomb blockade. The two sets of current- voltage characteristics are vertically shifted with respect to each otner 
by I/Ic = 0.3 for sake of clarity. 



6. Conclusions 

We have studied the current-voltage characteristics of mesoscopic Josephson junctions focusing on 
the modifications of the supercurrent as the junction parameters change. We have not discussed here 
the effect of the charging energy on quasiparticle tunneling relevant at higher applied voltages only 
pSfl . As the charging energy E c grows relative to the Josephson energy Ej , the supercurrent of the 
classical Josephson effect was shown to evolve gradually into a supercurrent peak caused by inco- 
herent Cooper pair tunneling. In the Coulomb blockade regime two types of structures may appear 
for small junctions with Ej <C E c embedded in a standard low impedance environment. The first 
structure, a peak at low voltages has recently been seen in experiments on lithographically fabricated 
junctions |ll], [l6| as well as break junctions p9| |. The second structure, peaks in the current-voltage 
characteristics due to resonances in the environmental impedance have also been seen in experiments 
[ p~5| , [l6] with a well-defined environment consisting of two transmission line segments allowing for a 
quantitative test of the theoretical predictions and good agreement was found. Although detailed 
experimental studies of the region between the phase diffusion and Coulomb blockade limits are 
absent, recent work Q indicates quantum effects in qualitative accord with the predictions made. 

We would like to thank M. H. Devoret, D. Esteve, and C. Urbina for stimulating discussions. One 
of us (HG) was supported by the Deutsche Akademischer Austauschdienst (DAAD), while support 
for the other author (GLI) during a stay at the Centre d'Etudes de Saclay was provided by the 
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